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presence of stray fields near the surface leads to an additional contribution to the repulsive 
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semi- infinite flux- line array. The presence of the superconductor/vacuum interface yields 
surface contributions to the compressional and tilt elastic constants. 
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1. Introduction 

The nature of the ordering of the magnetic flux array in the mixed state of high- 
temperature copper-oxide superconductors has received considerable experimental and 
theoretical attention in the last few years. It has been shown that fluctuation effects 
are important in these materials and can lead to a number of new phases or regimes, in- 
cluding entangled flux liquids, hexatic flux liquids and hexatic vortex glasses [f}^. Most 
experiments probe the properties of the flux array indirectly by measuring bulk properties 
of the superconductors, such as transport, magnetization or mechanical dissipation. At 
present direct measurements of the microscopic order of the magnetic flux array are mainly 
limited to decoration experiments at low fields These experiments aim to extract 

information on the vortex line configurations in the bulk of the material by imaging the 
pattern of the magnetic fiux lines as they emerge from the surface of the sample. It is 
clear that to interpret the experiments and assess whether one can indeed consider the 
surface patterns as representative of vortex line configurations in the bulk of the sample, 
one needs to quantitatively understand what are the relative effects of bulk versus sur- 
face interactions and disorder in determining the configuration of the vortex tips as they 
emerge at the surface. Almost thirty years ago Pearl showed that the interaction between 
fiux- line tips at a superconductor- vacuum interface decays as l/r±_ at large distances, with 
r±_ the distance beween fiux tips along the surface 0. In contrast, the interaction between 
fiux-line elements in bulk decays exponentially at large distances. For this reason Huse 
1^ recently questioned the assumption that is made implicitly in all experimental work 
that surface patterns are representative of flux lines configurations in the bulk. A proper 
interpretation of the fiux decoration experiments clearly requires a detailed knowledge of 



the modification of the properties of fiux-line arrays near the surface of the sample ||1C1|| . 
The presence of the surface modifies the pair interaction between fiux lines and changes 
the long wavelength properties of fiux-line liquids and lattices. 

The motivation of this paper is the desire to provide a framework for a quantitative 
analysis of the decoration experiments. The central result of the paper is a coarse-grained 
hydro dynamic free energy that describes the long wavelength properties of fiux-line liquids 
in a semi-infinite anisotropic superconductor occupying the half space ^ < 0. In a forth- 



coming publication ||T^ this free energy will be used as the starting point for evaluating the 



structural properties of the fiux-line tips as they emerge at the surface of a superconducting 
slab. 
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Explicit expressions for the liquid elastic constants, the compressional modulus and 
the tilt modulus, are obtained in terms of the superconductor parameters. The wave 
vector-dependent elastic constants can be written as the sum of bulk and surface contri- 
butions. The bulk parts are the well-known nonlocal elastic constants of a bulk flux-line 
liquid. In addition, the modification of the flux-line interaction near the surface and the 
magnetic energy from the stray fields in the region above the superconductor /vacuum in- 
terface lead to a surface contribution to both the compressional and the tilt moduli. These 
surface contributions vanish exponentially with the distance from the interface. They do. 



however, affect the surface properties and structure of the flux array [|Tl|. The compres- 
sional modulus of the flux-line liquid is enhanced in a surface layer near the interface by 
the stiffening of the repulsive interaction between the lines. At the interface the repulsive 
interaction between the flux tips decays as l/r±_ at large distances. This yields an additive 
surface contribution to the wave vector-dependent compressional modulus that diverges as 
l/qi_ at small wave vectors. In contrast, the surface contribution to the tilt modulus of 
the flux line liquid is always negative. This is because the magnetic fleld lines associated 
with each vortex fan out as the interface is approached. As a result, transverse magnetic 
field fiuctuations associated with tilting the flux lines become less costly in energy. The 
presence of the superconductor /vacuum interface produces a practically incompressible, 
but very flexible surface layer of flux-line liquid. 

Our flrst step to obtain the hydrodynamic free energy is the calculation of the mag- 
netic energy of an assembly of flux lines in a semi-inflnite superconductor sample in the 
London approximation. This calculation follows previous work by Brandt fl^ on vortex in- 
teractions near the surface of isotropic superconductors. Here we have in mind applications 
to the Cu02 superconductors and we consider the case where the interface is orthogonal 
to the c axis of the material. The magnetic fleld is also applied along the c axis, which is 
chosen as the z direction. The magnetic fleld energy is obtained in Section 2 by solving 
the London equation in the superconductor half space and Maxwell's equation in vacuum 
with the appropriate boundary conditions. The contribution to the magnetic energy from 
the superconductor half space is rewritten in a standard way as a pairwise additive inter- 
action between flux lines. The contribution from the vacuum can be naturally recast in 
the form of a surface modification to the pair interaction between the vortices. Similar 
results for the flux-line interaction in a semi-inflnite anisotropic superconductor were re- 
ported recently by Buisson et al. [0. No details were given in that paper. Since we think 
the derivation itself of the interaction is instructive and illuminating, we briefly sketch the 
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derivation here. In addition, the pair interaction is presented here in a form that, while 



equivalent to that of Ref. [jT^, is more transparent and more suited for approximations 
||llj| . Using this flux- line energy as the starting point in Sections 3 and 4 we obtain the 
coarse grained hydrodynamic free energy that describes the long wavelength properties of 
a semi-infinite flux array. The modification of the compressional and tilt elastic constants 
due to surface effects are identified and discussed. 



2. Magnetic energy of the flux-line array 

We consider a semi-infinite Cu02 superconductor sample that occupies the half space 
z < 0. The c axis is normal to the interface (i.e., along the z direction) and the sample 
is placed in a constant field H directed along the c axis, with Hd << H « Hc2- 
The local magnetic induction h(r) can be found by solving London's equation in the 
superconductor [z < 0) and Maxwell's equations in vacuum [z > 0) with appropriate 
boundary conditions at the surface = 0). The equations to be solved for an anisotropic 
uniaxial superconductor are: 



(2.1) 



(2.2) 



z<0: h + Vx(A-Vxh): 
V-h = 0, 

z>0: Vxh = 0, 
V ■ h = 0, 
with the boundary condition, 

z = : h continous 

(2.3) 

[z-(Vxh)],=o- =0. 

The fiux line is parametrized in terms of the coordinate z and Hi{z) = (ri(z)^ z) denotes 
the position of a fiux-line element at a heigth z below the planar superconductor/vacuum 
interface. The "line trajectories" ri{z) are assumed to be single- valued functions of z since 
"backtracking" involves a large energy cost for H > Hd. Below we will describe the 
tranverse fiuctuations of the fiux lines in terms of a three-dimensional line tangent vector 
Ti, defined as, 

T.iz) = ^^={Uz),z), (2.4) 
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or in terms of its component tj in the plane normal to the applied field, where, 

= (.5) 
The elements of the symmetric tensor A are 

= ^Ib^ap + (Ac - Xlb)CaCf3, (2.6) 

where Xab and Ac are the penetration lengths in the ab plane and along the c direction, 
respectively - c is a unit vector along the c axis. In the Ginzburg-Landau theory the 
anisotropy is accounted for by introducing an effective-mass tensor for the superconduct- 
ing electrons. For the high-Tg copper-oxides the effective- mass tensor is, to an excellent 
approximation, diagonal in the chosen coordinate system. Denoting by Mab and Mc the 
effective masses describing the interaction in the ab plane and in the c direction, respec- 
tively, the anisotropy ratio is defined as 7^ = Mc/Mab, and A^/A^^ = 7^- In the high-Tc 
materials 7^ >> 1. 

The boundary conditions are simply the continuity of the field at the boundary and 
the condition that there is no current normal to the interface. 

The solution of the equations can be obtained directly in terms of the partial Fourier 
tranform of the local induction, 

h(q±,^) = j dv^e-'^^-^^\i{v), (2.7) 

where r = (r_L, 2;), with the result, 

h(q^,^) = ^y r dzie-^q^-'-^("^) ^~'^^ e"^"e-^^", (2.8) 
Kb i J -00 g± + a 

for 2; > 0, i.e., in the vacuum half-space, and 



h(qi,^) = -^$^ f dzie-^^-'-^(^^)|[z + q^(q^-t,(^i)) 



g-Q:|Z-Zi| 



a 



p -, ga(2 + 2l) 2e"(2 + 2i) 

+ z - q_L(q_L • ti(2;i)) [zq-i + zq± a] — ■ — - (2-9) 

L ^ ' \ a aiq 1 -|- o;) 



+ (z X q^) (z X q±) •ti(2;i) 



a{q±_ -\- a) 

g— ac|2— zil _ gac(2:+2:i) 



for 2; < 0, where 



a = a{qi_) = ^ql + 1/A^, , (2.10) 



and 

«c = «(79±) = \[i^^+VK, ■ (2-11) 

Alternatively, one can find the local induction by using the method of images, as was 
done by Brandt for an isotropic superconductor [|12|. If the line is straight, the field of 
a single flux line is identical to that evaluated many years ago by Pearl The field of 
many vortices is just the linear superposition of the field of a single vortex. 

The total magnetic energy is given by, 

U = U^ + Us= [ dr^ + ! dr^ [h^ + {V x h) ■ A ■ {V x h)l , (2.12) 

Jz>o Stt J^<o Stt 

where the first term on the right hand side is the field energy from the vacuum half 
space and the second term is the energy from the fields and the supercurrents in the 
superconductor. The vacuum contribution Uy is immediately evaluated by substituting 
Eq. ( p.8|) in the first term of Eq. ( |2.12|) , with the result. 



The contribution Ug from the half space occupied by the superconductor can be evaluated 
either by direct substitution of Eq. (|2.9| ) or by first performing an integration by parts to 
obtain. 



Us 

8n 



-[ rfr h- fh + V X (A- V X h)l + ^ /"rfr^ z- [h X (V X h)l (2.14) 

In J^<o L i Stt J ^^-u 



By substituting from the London equation (|2.1|) in the first term on the right hand side of 
Eq. ( p.l4|) , one then obtains, 

^^ = ^E / dziT,{zi)-h{r,{zi),zi) + ^ y"(irxz-[hx(Vxh)]^^Q_, (2.15) 

where h.(ri{zi), Zi) is the local field due to all the fiux lines, as given in Eq. (|2.9|), evaluated 
at the location of the z-th line. After some manipulation of the surface contribution in Eq. 
( p.l5|) , the total energy can be written as an integral over the superconductor half space 



of a pair interaction between flux lines, 



U 



SttA 



X2. 



dzi 



(iq 



(271) 



^-=-«qx-[ri(2i)-rj(22)] 



-Ol\zi-Z2\ 



2a 



■(1 + ti-tj 



+ 



+ 



-ac\zi—Z2\ 

2a, 

^a{zi+Z2) 

2a 
z X q±) 



-a zi— 22 



2a 



1 - 1, 



(z X qx) ■ t 
2 



(z X q_L) ■ tj 



ti (z X q^) ■ tj 



Q — ac{zi+Z2) ga(2:i+22)-| 



2ar 



2a 



(2.16) 

where = ti{zi) and tj = tj{z2). The energy of Eq. ( |2.16D includes the self-energy 
of the vortices (terms with i = j in the sum). The small distances divergence of the 
repulsive self-energy needs to be truncated by introducing a cutoff to account for the 



finite size of the vortex core. As discussed by Sudbo and Brandt the proper cutoff is 
anisotropic. The terms in Eq. ( p.l6| ) that contain exponentials in 1^1 — ^21 are identical to 
those obtained for the pair interaction between flux lines in bulk. In fact if one assumes 
that the superconductor extends over all space in the z direction, it is easily shown that 
the part of the interaction containing exponentials in \zi — Z2\ is simply identical to that 
obtained elsewhere for a bulk superconductor |[1^,|16[ . The terms containing exponentials 
in {zi + Z2) arise from the fluctuations in the local induction due to the presence of the 
superconductor /vacuum interface. The corresponding contribution to the energy density 
is nonzero only within a layer of depth ^ Xat near the surface. The total energy ( |2.1(j| ) 
of the flux array also includes the magnetic energy from the stray fields generated by the 
vortex ends outside the sample. 

It is useful for the following to rewrite the interaction energy (|2.16|) in a form where 
the transverse (to qj^) components of the tangent vectors are eliminated in terms of the 
magnitude of the vectors and their longitudinal components. This can be done by using 
[(z X q_|_) ■ ti] [(z X q_|_) ■ tj] = ■ tj — (qj_ ■ t^) (q^ ■ t^). The terms containing only the 
longitudinal components of t^ are then integrated by parts, with the result. 



U 



8nX 



ab ij 



dzi 



dzo 



(iq_L 



-ici±-[ri{zi)-rj{z2)] 



(27r)2 

Zl-Z2\ _|_ gQ(2l+22) 



2a 

30(21+22) 



2 



^—ac\zi-Z2\ _|_ gac(2l+22) 



^ (2.17) 



-ac\zi-Z2\ 



+ 



yab(l^oi{q± + a) 



-(- tj ■ ti 



2ac 

^0^(21+22) 



2ar 



This form of the energy is equivalent to that given in Eq. ( |2.16|) above and is a convenient 
starting point for the derivation of the elastic energy of a flux array described below. 

To gain some physical insight on the effect that the presence of a superconduc- 
tor/vacuum interface has on the flux-line interaction, it is useful to consider the interaction 
for the case of straight flux lines. In this case the two-dimensional position vectors do 
not depend on z and the tangent vectors simply vanish. The integration over z\ and 
can be carried out. The interaction energy of an array of rigid flux lines is then given by, 



-\Y.\ tSt^^""^'""^^^ V + ^^^^^ 



?_La(?± + ol) 



(2.18) 



where L is the size of the superconductor sample in the z direction (L >> Agb) [pTTf. 
There are two contributions to the energy of Eq. ( |2.18| ). The first term is a bulk energy 
proportional to the size of the system in the z direction. The pair potential Vs((?±) is 
given by 

and it is the Fourier transform of the usual pair interaction per unit length between straight 
flux lines in bulk. The second term is a surface energy corresponding to the magnetic energy 
of the stray fields near the interface. It can be interpreted as a pair interaction between 
flux-line tips at the superconductor /vacuum interface, with 

VsiqA = — ^ , (2.20) 

In the long wavelength limit the surface interaction becomes 

Vsi.,.) « (2.21) 

Inverting the Fourier transform, one finds that the pair interaction between flux-line tips 
decays as Vsiri) ~ </'o/4^^^_l at large distances (r^ >> Aat), as obtained many years 
ago by Pearl 0. As pointed out by Huse 0, this is easily understood because each flux 
line spills a flux quantum (^o iiifo the vacuum half-space when exiting the superconductor's 
surface. The interaction between flux tips is therefore the interaction between two magnetic 
monopoles of "charge" (/)o/27r. 
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3. Elastic Energy 

In this section we calculate the continuum elastic energy associated with long- 
wavelength deformations of the semi-infinite flux array. We confine ourselves to applied 
fields in the range Hd << H << ifc2, corresponding to Xab > d >> ^a&- Here d = ^/no 
is the average intervortex spacing in the xy plane, with tiq = (^q/Bq the equilibrium 
areal density of flux lines and Bq the equilibrium induction in the superconductor. For 
H » Hci, Bq ~ H. For the flelds of interest the vortex cores do not overlap and the 
spatial variations in the order parameter outside the core can be neglected. The energy of 
an array of London vortices that was obtained in Section 2 is then the appropriate starting 
point for obtaining the continuum elastic energy of the flux array. 

In order to calculate the energy associated with elastic distortions of the flux array, 
we follow Ref. and write the two-dimensional position vector of the flux lines as. 



ri{z) = Tieq + Ui(z), (3.1) 

where r^eq are the equilibrium positions and Ui{z) two-dimensional displacement vectors 
in the plane normal to the applied fleld. Strictly speaking the equilibrium positions in Eq. 

should be the equilibrium positions of the flux lines in a semi-inflnite superconductor 
in the presence of a surface. In general the equilibrium solution of the semi-inflnite problem 
will differ in a surface layer from the usual Abrikosov solution of the bulk problem |T^] . On 



the other hand, we are interested here in evaluating the long wavelength elastic energy in a 
regime where d < Xab and the magnetic flelds of the vortices overlap. We will then neglect 
below all corrections to the elastic energy due to the discreteness and the speciflc structure 
of the flux-line lattice and replace all lattice sums by integrals. It is then consistent to 
also neglect the deviations of the equilibrium flux-line positions near the surface from a 
regular Abrikosov lattice. As a result of this continuum approximation the elastic energy 
obtained below contains compressional and tilt elastic constants, but no shear modulus. 
To evaluate the shear energy one needs to carry out a more microscopic calculation that 
incorporates explicitly the discreteness of the flux lattice. Since our ultimate interest is in 
obtaining the hydrodynamic free energy of a semi-inflnite flux-line liquid, the calculation 
of the shear modulus is beyond the scope of this paper. 

We then assume that the equilibrium positions r^eq are everywhere those of a regular 
triangular Abrikosov lattice and expand the energy of Eq. ( p.l7|) for small displacements 



from equilibrium, retaining terms up to quadratic in the displacements. The details are 
sketched in Appendix A. The terms linear in vanish and one obtains, 



[/ = C/eq + (3.2) 



with 



^U = ^Y.f r ^^2^ J]|[(q± + Q).(q± + Q)/3Gi(|q^ + Q|;zi,Z2) 

- QaQf3Gi{Q;zi,z2) Wa(q±,^i)M/3(-q±,22) '^^•^^ 

+ G'2(q±; Zi,Z2)dz^u{q±, zi) ■ d:,2u{-q±, Z2) 
where 

Ci g±;2i,22 =T^i ^ + ^ T 

Q:(zi+Z2) 



+ 



A„feg±a(g± + a) 
and 



(3.4) 



/,2 p—aa\zi—Z2\_^oia{zi+Z2) 



G2{q±; zi, Z2) = . (3.5) 

^^Kb 2ae 

Here Q are the reciprocal vectors of the triangular Abrikosov lattice and A is the area of 
the system in the plane normal to the applied field. Also, u(q_L,2;) is the lattice Fourier 
transform of the displacement, as defined in Appendix A. 

Since we are only interested in the continuum limit here, we will not attempt to 
perform explicitly the sum over the reciprocal lattice vectors. To proceed, we separate out 
in Eq. ( |3.3| ) the Q = term in the sum over Q. This term gives the collective contribution 
to the energy that dominates for long-wavelength elastic deformations. In the remainder 
of the elastic energy, containing XIq^c neglect qj^ compared to Q. This is because 
IQI > ksz, for Q 7^ 0, where kBz = \/2Tcno is the size of the first Brillouin zone, and we 
are interested in deformations of the flux array with q± << ksz- The elastic energy can 
then be rewritten in a form that explicitly identifies a compressional and a tilt energy, 

5f/ = — ^/ dzi dz2<B{q±;zi,Z2) ci±-u{(i±,zi) ci±-u{-q±,Z2) 

^(<?±; zi, Z2)d;,^u{ci±, zi) ■ d:,^u{-ci±, Z2) 



where B{q±; zi, Z2) and }C{q±; zi, Z2) are the compressional and tilt elastic constants per 
unit length, given by, 



2a V 



^—ac\zi-Z2\ \ ^ac{z\+Z2) q,2 ^01(21+22) 
+ ^ ^ + 



(3.7) 



2ac gi <?±A^^a(g± + a) 

and 

^2 g-ac|2l-22| _ gQc(2l+22) 

^{qr, zi, Z2) = ICq^zi, Z2) + . (3.8) 

In this continuum approximation the Q 7^ part of the sum in Eq. (^]^) vanishes when the 
displacements do not depend on z, i.e., the lines are rigid. Therefore it does not contribute 
to the compressional modulus, but only to the tilt energy. The corresponding contribution 
is denoted by /Co(^i5 ^2) and is independent of qj^ because in this term we have neglected 
q_L compared to Q > ksz ■ It is given by 



M^i,^2)=^^):^| 

ab o^n 



Bn(t)n 1 ^ f e-"(^)l^i-^2| _ ea(Q)(2i+z2 



+ 



2A>3(Q) 

g-ac(Q)i2l-22l _ gac(Q)(2l+22) 



2ac(Q) (3.9) 

(Q)(2l+22) 



e 



A>^(Q)[Q + «(Q)] 

- e(zi - Z2)e^2"(^) - e(z2 - ^i)e^i"(^) 

In the limit where the equilibrium areal density of vortices goes to zero, /Co(-Si, -22) reduces 
to the tilt contant of a single flux line. 

The elastic constants of a semi-infinite superconductor usual, nonlocal in the 

xy plane and contain two types of nonlocal effects in the z direction. As in a bulk sample, 
the elastic constants depend on the vertical distance \zi — ^2! between any two small 
volumes of the elastic fiux array. This nonlocality refiects directly the range of the repulsive 
interaction and it occurs on the scale of the penetration lengths. In addition, the elastic 
constants of a semi-infinite superconductor depend on the distance of each deformed flux 
volume from the superconductor/vacuum surface. These surface effects yield the terms 
that depend exponentially on the distance of the deformed flux volume from the surface - 
the exponentials in (-Zi +-^2) or zi and Z2 in Eqs. ( |3.7| ) - ( |3.9| ) . They are important within a 
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surface layer of thickness determined by the penetration lengths. In Appendix B we display 
the expression for the elastic energy obtained by replacing the integrals over zi and Z2 by 
wavevector sums in the corresponding Fourier space. This expression is instructive because 
it contains wave vector-dependent elastic constants that naturally separate into the sum 
of the well-known wave vector-dependent elastic constants of a bulk flux-line array and 
surface contributions. 

In order to gain some physical insight on the surface modification of the elastic con- 
stants of the flux array, it is useful to consider the elastic energy corresponding to two 
specific deformations: an isotropic compression and a uniform tilt of the flux array. 



Isotropic Compression 

Consider a pure isotropic compression of the flux array, where ■ u 7^ 0, but u is 
independent of z. The corresponding elastic energy dUcomp is immediately obtained from 

m) as, 

SUcomp = -^'^ dzi dz2 B{q_L;zi,Z2) \ci± ■u{ci_l)\'^. (3.10) 



The z integrations can be carried out in Eq. (|3.10|) , with the result, 

SU^omp = Bsiq^) + B^iq^)] ■ u(q^)|2, (3.11) 

where 

S3(?±) =CL(g±,g. = 0) = ^ \ (3.12) 
is the compressional modulus of a bulk flux-line array, as given in Eq. (B.5) and obtained 



before by other authours (see, for instance, Ref. [|T^,|T^), evaluated at q^ = 0, and 

B2[q±) = — 72 — 7 — ^ (3.13) 

is the compressional modulus of the two-dimensional array of flux-line tips at the super- 
conductor's surface. In the long wavelength limit, i.e., for q±Xab << 1, B2{q±) reduces to 
the two-dimensional bulk modulus of an array of monopoles, interacting via a 1/r^ poten- 
tial, B2{q±) — Bs{q±) /q± ^ ^o/(47r^q'±) @. As discussed earlier the repulsive interaction 
between flux lines becomes stronger and long ranged near the superconductor/vacuum 
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interface. This yields a corresponding increase of the compressional energy of the flux 
array. 



Uniform Tilt 

We now consider the energy corresponding to a uniform tilt of the flux array, that is 
a deformation with (z x q^) ■ dzu{q±, z) = 9{q±), but ■ u = 0. The corresponding tilt 
energy is 



tilt 



-y 

2A ^ 



dzi dz2 IC{q±;zi,Z2)\9{q±] 



or, carrying out the integrations over the z variables. 



(3.14) 



(3.15) 



where K3{q±) = CM{q±,qz = 0) is the tilt coefficient of a bulk flux array, given in Eq. 
(B.6) and obtained before [^|TB|,|T^, evaluated at q^ = 0, 

1 



4^1 + ?iA2,7^ 



+ C44(?z = 0), 



with 



C44(?z = 0) = nof— ^ 
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^ ln(7«) - A In (1 +7'A^,4^) + 



2Y 



2(1 + K,klz 



(3.16) 



(3.17) 



and K2{qi_) is the surface contribution to the tilt constant, 
K2{qi.) 



BoXab 1 

47r (l + ?iA2,72)3/2 

+ noXab(-r-T— 

\4nXab 



1 



2 

+ 77- 



1 



1 



hW^ + >^'iA 3 (1 + A2,fc|^)3/2 
1 



(3.18) 



+ 



1 



— tan 



ab'^BZ 
1 



Xabksz + + X'i^kl 



•'BZ ^^abr^BZ -r V ^^ab'^BZ ' 

The surface contribution to the tilt energy is always negative. This can be understood 
physically because the magnetic field lines associated with each vortex spread out as the 
interface is approached. The magnetic field fluctuations associated with each line are 
appreciable in an area that becomes very large near the interface. Consequently the energy 
cost for tilting the lines or inducing transverse line fluctuations decreases. 
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4. Flux Liquid Free Energy 

In the range of applied fields of interest here, Hd « H « Hc2, the properties of a 
flux-line liquid on length scales large compared to the intervortex spacing can be described 
in terms of two hydrodynamic fields, a microscopic areal density of vortices, 

n„,ic(rx,^) =5^5^'^(r±-r,(^)), (4.1) 

i 

and a microscopic "tangent" field in the plane perpendicular to the applied field, 

tmic{v^,z) = 5^^\v^-v,{z)). (4.2) 

i 

From these microscopic densities one constructs coarse-grained density fields n(rx, z) and 
t{r^,z) by averaging ( [4.1|) and (|4.2|) over a hydrodynamic volume centered at rj^. As in 
the Landau's theory of phase transitions the long wavelength properties of the fiux-line 
assembly can be described in terms of a coarse-grained free energy retaining only terms 
quadratic in the deviations 5n(rj^,2;) = n(rx,z) — no and t{vj^,z) of the hydrodynamic 
fields from their equilibrium values. Because we are dealing with magnetic fiux lines that 
cannot start or stop inside the medium, the density and tangent fields are not independent, 
but satisfy a "continuity" equation in the time-like variable z, 

a^n 4- ■ t = 0. (4.3) 

This condition refiects the requirement of no magnetic monopoles. It can be implemented 
by introducing a vector potential or two-component "displacement field", u(rj^, z), with, 

dn = -noV_L • u, (4.4) 

and 

du 



It was shown earlier [ffUl that when the density and tangent fields are expressed in terms of 



the displacement vector u, the hydrodynamic free energy of a fiux-line liquid differs from 
the continuum elastic free energy of the Abrikosov fiux-line lattice only for the absence 
in the former of a shear modulus. In other words it was found that in a bulk sample the 
hydrodynamic fiux-line liquid free energy, when expressed in terms of the vector potential 
u, is identical to the continuum elastic energy of a fiux-line lattice, provided all effects due 
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to the discreteness of the lattice - in particular the shear modulus - are neglected in the 
latter. Similarly, one can show that the compressional and tilt moduli of a semi-infinite 
fiux-line liquid are identical to those obtained in Section 3. The hydrodynamic free energy 
of a fiux-line liquid is then given by 

j^r^ dzi dz2<B{q_L;zi,Z2) 5ny{ci_L,zi) 5ny{-ci_L,Z2) 
Zi,Z2)t{q±, zi) ■ t(-q^, Z2) k 



where B{q±; zi, Z2) and JC{q±; zi, Z2) are the compressional and tilt elastic constants per 
unit length, given in Eqs. ( ^77| ) and ( p.8|) , respectively. The statistical averages over ( ^?6| ) 
need to be evaluated with the constraint ( |4.3| ) . If the constraint is incorporated in the free 
energy by expressing density and tangent fields in terms of the vector potential u according 
to Eqs. ( |4.4| ) and ( [4.5| ), the flux-liquid free energy and the elastic energy of the flux-line 
lattice only differ for the absence in the former of the shear modulus. 

The flux-line liquid free energy given in Eq. ( [4.6| ) can also be obtained directly from 
the flux-line energy of Eq. ( p.l6|) by the coarse-graining procedure described for instance 



in Ref. |]T^, without making any references to an equilibrium lattice of flux lines. Care 
has to be taken in dealing with the self-energy term, corresponding to the i = j part of 
the pair interaction. 

In a forthcoming paper we will propose an approximate form of the hydrodynamic 
free energy ( [4.6|) obtained by assuming that the most importance source of spatial inho- 
mogenieties in the z direction is the presence of the surface itself and by neglecting all 
other nonlocalities in the z direction. This approximate hydrodynamic free energy will be 
used there to analyze the interplay of bulk and surface forces in determining the structure 
of the flux-line tips as they emerge from the superconductor sample. 
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12330. It is a pleasure to thank D.R. Nelson for stimulating discussions and the Diparti- 
mento di Fisica, Universita di Roma "Tor Vergata" , for hospitality during the completion 
of this manuscript. 
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Appendix A. Derivation of the Elastic Energy 

When the energy given in Eq. ( |2.16| ) is expanded for small displacements of the flux 
lines from their equilibrium positions, the first nonvanishing correction to the equilibrium 
energy is quadratic in the displacements and it is given by, 



X 



+ G2{q±; zi, Z2)d^^Ui{zi) ■ d^^Uj{z2)^, 



where the elastic kernels Gi{q±] zi, Z2) and G2{q±] zi, Z2) Hve been given in Eqs. ( p.4|) and 
( p.5|) . It is convenient to expand the lattice displacements in Fourier series according 
to, 

u,(^) = - u(k,^)e^'^-'--, (A.2) 



an 

^ k€SZ 



where ao = l/n-o is the area of the primitive unit cell and the sum is restricted to the 
wavevectors of the first Brillouin zone, as indicated. Conversely, the Fourier coefficients 
u(k, z) are given by 

u(k, ^) = — V u,(z)e-^''-'"-'^. (A.3) 
no ^ 

The normalization of the Fourier expansion ( |A.2|) has been chosen in such a way that 
the lattice Fourier amplitudes u(k, z) as defined in Eq. ( [A3|) for kGSZ are also the 
two-dimensional continuum Fourier transform of a coarse-grained density field, 

u(r^,^) = — Vu,(;2)5(2)(ri -r,eq), (A.4) 
no ^ 

according to. 



u(k,z) = / rfrie-^''-'"^u(ri,2). (A.5) 



The (5-function in Eq. (|A.4| ) is really a smeared-out two-dmensional 5-function with a finite 
spatial extent ksz- After inserting the Fourier expansion ( |A.2| ) in Eq. ( |A.1| ) , the lattice 
sums can be carried out, using. 



Q 
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where Q are the reciprocal lattice vectors of the Abrikosov lattice, and 



(A.7) 



keBZ 



Finally, making use of the periodicity of the displacements in reciprocal space, u(k + Q) 
u(k), one obtains Eq. (|3.3| ) . 



Appendix B. Wavevector-Dependent Elastic Constants 



It is useful to rewrite the elastic energy of Eq. ( |3.6| ) by taking the Fourier transform 
of the displacement vectors with respect to z, 



(B.l) 



u(q±,gz)= / dz e "^^'^u(q_L, 2;), 
J —00 

for Im[qz) > 0. The elastic energy can then be written as 
^ ^ dQz dq' { 

Qz, qi)qzq'M<l±, Qz) ■ u(-q±, q'^[ 



The wave vector-dependent elastic constants are given by. 



B{q±] qz, q'z) =2vr5(Q, + q'^) CL{q±, q^) + <^ ( 



1-^ 



1 



1(lz 



+ 



+ 



q^^ J a{a - iqz) + c^^ 
1 1 



(f^ioLc- iqz) o? (^z + g_L) {a - iqz){a - iq'^ 



and 



1 



^Q±; Qz, q'z) =271(5(9^ + q^) C44(?±, - -^—^2 ■ 2^ — 12 

'^^Kb «c - iQz ai + q'/ 

1 1 



+ 
+ 
+ 



47rA?^ 2A 
1 



1 



ac{Q)c - iqz aliQ) + q 
Q 



1 



1 



[Q + am La(Q) - t{qz + q'z) Va(Q) - % 
1 \ 1 



a{Q)~iq'z^ {a{Q) -iqz){a{Q) -iq'^ 
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fB.2) 



(B.3) 



(B.4) 



In this form the elastic constants are exphcitly given by the sum of bulk and surface 
contributions. The bulk contributions are the usual ones, given by, 



CLiq±,qz) - 



(B.5) 



with q"^ = q']_ + q^, and 



C44(<?±,?z) 



+ C44(<?z)- 



(B.6) 



In Eq. ( p.6| ), 044(^2) is the contribution to the tilt coefficient arising from the large 
wavevector (Q 7^ 0) part of the lattice sum. The sum over the reciprocal lattice vectors is 
evaluated in the continuum limit by replacing it by an integral with appropriate cutoffs, 
according to, 

1 . — . /' rin 

(B.7) 



-y 



kBZ<Q<\l^o.b 



f27r)2- 



For K = Xab/iab >> 1 oue obtains. 



C44(?2) =no (4^-^) I ln(7K) - ^ In (l + j^Xl^^klz + >^lbql) 



+ 



2 ^2 



In (1 + Xl^,klz + Xl^ql) - In (1 + Xl^klz) 



Aside from numerical differences due to the details of the short-length scale cutoff, the 
expression for C44 obtained here agrees with that given by D.S. Fisher and, as discussed 
there, corrects an earlier result of Brandt and Sudbo . Finally, in the limit where the 
density no of flux lines vanishes - and therefore ksz — > in Eq. ( [B.8| ) -, C4^4^{qz) / reduces 
to the single-line tilt coefficient, ei{qz), given by 



ei(<?z) = lim 044(^2) /no 

no^O 



= ( \^ 

For = 0, one obtains 



1 

J2 



ln(7K) 



27 



^ ln(l + A2,g2) + _J_in(i + A2,Q2^ 

2 



(B.9) 



2 r 1 1 



:b.io) 
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